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>- ■ Metal line ratios in a sample of 13 quasar spectra obtained with the HIRES 

spectrograph on the Keckl telescope have been analyzed to characterize the evo- 
lution of the metagalactic ionzing flux near a redshift of 3. The evolution of 
Si IV/C IV has been determined using three different techniques: using total col- 
umn densities of absorption line complexes, as in Songaila & Cowie (1996); using 
the column densities of individual Voigt profile components within complexes; 

Q_i' and using direct optical depth ratios. All three methods show that Si IV/C IV 

changes abruptly at z ~ 3, requiring a jump in value of about a factor of 3.4, 
and indicating a significant change in the ionizing spectrum that occurs rapidly 
between z = 2.9 and z — 3, just above the redshift at which Reimers et al. (1997) 
detected patchy He II Lya absorption. At lower redshifts, the ionization balance 
^ ■ is consistent with a pure power law ionizing spectrum but at higher redshifts the 

spectrum must be very soft, with a large break at the He + edge. An optical 
depth ratio technique is used to measure the abundances of ions whose transi- 



tions lie within the forest and C III, Si III and O VI are detected in this way. The 
presence of a significant amount of O VI at z > 3 suggests either a considerable 
volume of He III bubbles embedded in the more general region where the ioniz- 
ing flux is heavily broken, or the addition of collisional ionization to the simple 
photoionization models. 

Subject headings: early universe — intergalactic medium — quasars: absorption 
lines — galaxies: formation 
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1. Introduction 

Whereas we know from the absence of any significant Gunn-Peterson effect even in 
the highest redshift quasars that hydrogen reionization of the intergalactic gas must have 
taken place at z > 5, we have much less information about the period at which the bulk 
of singly ionized helium converted to doubly ionized helium. Since late He + ionization may 
significantly change the temperature of the intergalactic gas it is critical to understand this 
heating if we are to correctly model the growth of structure in the IGM, and determine 
the mapping of the baryon density to observable quantities such as observations of the 
neutral hydrogen Lya forest and the He + Lya opacity. Phenomenological modelling of this 
event depends critically on the softness of the composite spectrum of the ionizing sources 
(e.g. Miralda-Escude & Rees 1993; Madau & Meiksin 1994), amplified by the subsequent 
radiative transfer, and such models cannot be considered reliable in predicting the high 
energy (E > 54 eV) metagalactic ionizing spectrum above the He + ionization edge, which 
determines the fraction of singly ionized helium. 

Our most direct information on the He + opacity is through observations of quasars 
whose spectra extend to the He + Lya wavelength. Despite the extreme difficulty of these 
measurements, successful observations of the He + Lya absorption have been made toward 
z > 2.8 quasars with HST (Jakobsen et al. 1994; Hogan et al. 1997; Reimers et al. 1997) 
and of the z = 2.72 quasar HS1700+6416 with HUT (Davidsen et al. 1996; Zheng et 
al. 1998). The He + Lya opacity shows a marked decrease from a value of r = 3.21^ in 
Q0302-003 at z = 3.29 to r = 1.0 ± 0.07 in HS1700+6416. More remarkably, the Reimers 
et al. observation of the intermediate redshift quasar HE2347— 4342 at z — 2.89 shows both 
'troughs' and 'voids' in the He II Lya observations, suggesting that at z ~ 2.8 we are seeing 
fully ionized He III bubbles interspersed among as yet unionized He + regions and that it is 
at this point that the porosity of He III regions is approaching unity. 

The recent discovery that the bulk of Lya forest absorption with iV(H I) > 3 x 10 14 cm -2 
contains associated metals (Songaila & Cowie 1996, hereafter SC) gives us an alternative 
approach to the problem since ionization balance in these forest metals provides a diagnostic 
of the shape of the metagalactic flux in the neighborhood of the He + edge. As was first 
noted in Songaila et al. (1995), the value of Si IV/C IV in high ionization systems is 
critically dependent on the He + ionization edge break strength. This has subsequently been 
investigated in more detail by SC, Savaglio et al. (1997) and Giroux & Shull (1997), among 
others. As will be discussed further here, the bulk of the forest metal line systems at the 
currently observed redshifts (z ~ 2 — 4) are high ionization (C II/C IV <C 0.1) and so, unless 
there is a strong break at the edge, they will have Si IV/C IV 0.1 even for the higher Si/C 
abundances characteristic of low metallicity systems. Therefore, once the He + is fully ionized 
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and the IGM becomes relatively transparent to the integrated quasar spectrum (e.g. in the 
metagalactic spectrum of Haardt & Madau 1996), the observed Si IV/C IV values should fall 
in this low range. SC and Savaglio et al. (1997) have shown that, whereas this is generally 
true at z < 3, much higher Si IV/C IV values are regularly seen at z > 3, suggesting a 
significant change above this redshift, which would be consistent with the interpretation of 
the Reimers et al. (1997) observations as showing the redshift at which He III bubbles begin 
to overlap. Boksenberg (1998) has recently questioned this result, based on an analysis of the 
the redshift evolution of the ion ratios in the separate Voigt profile components in complex 
systems rather than in the integrated column densities of the complexes. However, his 
analysis is based on rather a small sample of systems without clear selection criteria. In this 
paper, I shall use the largest sample to date to demonstrate unambiguously that, irrespective 
of the method of analysis, there is indeed a rapid jump in the value of SI IV/C IV at a redshift 
just below 3, and that the ionization stages in the metals are consistent with this being, for 
most systems, the point at which they change from being ionized by a metalgalactic spectrum 
that is, at z > 3, heavily broken above 54 eV to one that is only mildly broken at the lower 
redshifts. 

The sample and the data reduction are described in §2 and the reader who is primarily 
interested in the results could safely skip this section. The evolution of the C II/C IV 
and Si IV/C IV values with redshift is described in §3 where I show that the presence of 
a jump in Si IV/C IV values at z just under 3 is highly significant and does not depend 
on the method of analysis, whether by total column densities of complex, by the column 
densities of individual Voigt components, or by directly analyzing the distribution of optical 
depth ratios. In §4 I consider the overall ionization balance including intermediate ionization 
stages such as C III and Si III and high ions such as N V and O VI using, for those lines 
that lie primarily within the forest, the optical depth distributions of the ensembles of such 
lines, which provides a new and robust technique for determining their properties. The 
overall ionization balance is broadly consistent with unbroken power law photoionization at 
z < 3; however, the observation of significant amounts of O VI at z > 3 requires either the 
presence of a considerable volume of He III bubbles permeating regions where the ionizing 
flux is heavily broken, or the addition of collisional ionization to the simple photoionization 
models. Finally, the conclusions are briefly summarised in §5. 



2. Data 



The sample was drawn from observations of 13 quasars with complete or near-complete 
coverage between the quasars' Lyman alpha and C IV emission lines. All observations (sum- 
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marised in Table []]) were made with the HIRES spectrograph on the Keckl 10m telescope, 
using a 1.14" x 7" slit, which gives a resolution of R = 37, 000, and with total exposure times 
from 3.3 hours to just over 10 hours. The approximate S/N per resolution element for each 
quasar at wavelengths between Lya and C IV emission is summarised in Table []] and ranges 
from 60 (Q0256-000) to 330 (Q1422+231). 

Individual exposures were generally set to 40 minutes (providing a reasonable compro- 
mise between minimizing cosmic ray hits and read noise effects) and with the 2048 Tektronix 
CCD read out with a binning of two in the spatial direction and one in the spectral direction, 
resulting in a spatial pixel of 0.38" and AA/A = 7.3 x 10~ 6 ; this substantially oversamples 
the resolution element. For the higher redshift quasars several settings of the HIRES cross- 
disperser were used to provide the full wavelength coverage. All but three of the quasars 
have full coverage of the wavelength region of interest (Table 1, column 4) and the remain- 
ing three have small interorder gaps. In general, the spectra extend from approximately 
3800 A, below which the HIRES sensitivity begins to fall rapidly, up to the wavelength of 
the quasar's C IV emission. A white dwarf star and the thorium-argon calibration lamp were 
also observed either immediately before or immediately after each sequence of exposures on 
a quasar. 

A relatively straightforward extraction technique was used since optimization does not 
provide substantial gains for these spatially undersampled and signal-dominated spectra. 
Starting from approximate order and wavelength solutions, the position and shape of each 
order on the CCD was determined from the observations of the white dwarf. Only orders 
lying fully on the CCD were extracted. A small number of known bad points and columns 
were cleaned by interpolation, and a number of regions, such as the central blemish on the 
chip, flagged as unusable. The first pass cosmic ray rejection was made by forming a residual 
image with a 5-point two-dimensional median subtracted and then subtracting an 11-point 
row by row median. Any points in this differenced image that exceeded expected values 
by more than 4 a were interpolated from neighboring points along the column (spatial) 
direction. 

The spectral image was distorted by bilinear interpolation to make the echelle order 
parallel to the rows. The sky was then determined from the median of pixels on the same 
column lying between 2" and 3.5" on either side of the maximum signal, and subtracted 
from the spectral image, and the spectrum was extracted in a ±1" window. The spectra for 
all the sub-exposures were then combined with an optimal S/N weighting, and any residual 
cosmic ray structure was removed by identifying any strongly deviant (> 5 a) points in the 
individual spectrum and flagging and removing these from a second and final summation. 
The blaze function was determined from a normalised fifth-order polynomial fit to the white 
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dwarf spectrum, iteratively repeated to exclude any absorption lines, and divided into the 
spectrum. The wavelength and spectrum for each order were then individually stored. 

The Th-Ar calibration spectrum was extracted in the same fashion as the object spec- 
trum but without any sky subtraction step. For each order, each of the stronger unblended 
calibration lines was fitted with a gaussian to determine its centroid, and a third-order poly- 
nomial fit was then applied to determine the wavelength solution. From 10 to 30 calibration 
lines were used for each order and the r.m.s. residual was normally less than 2 x 10~ 3 A, so 
that velocities should be accurate to better than 0.1 km s -1 . Finally, the wavelength scale 
was converted to vacuum heliocentric, which is used throughout the present paper. 

A typical extracted order longward of the forest is shown in Figure [la| which shows 
order 63 from a 2-hr exposure of Q1159+123, obtained as three 40-minute exposures in 1997 
April. The absorption lines in the center of the order are a C IV doublet, and the rise in flux 
to shorter wavelengths reflects the onset of the quasar's Lya and N V emission. Figure |Tb" 



shows an order in the forest blueward of the quasar's Lya emission, where the accuracy of 
the sky subtraction may be gauged from the precision with which the bases of the saturated 
lines approach zero. 

A final combined spectrum was made from all the observations taken at all echelle 
settings, with spectra covering each order added in a S/N optimized fashion to form a 
combined spectrum of that order. The final 8-hr spectrum of order 63 in Q1159+123 is 



shown in Figure g for comparison with Figure [La]. An associated effective exposure time 
vector, normalised by the S/N, was also maintained at each point; this can be combined 
with the counts to determine the noise at each point in the spectrum. For some purposes it 
is useful to have a combined spectrum and this was formed by interpolating all the orders to 
a uniform (0.05 A) grid which provides Nyquist sampling throughout the wavelength range. 
A typical spectrum (of Q1159+123) is shown as Figure |3|. This spectrum has a Lyman 
limit system just longward of 4150 A and stretches through to C IV emission just below 
7000 A, near the atmospheric A band, which can be clearly seen in the spectrum. Most of 
the remaining features between Lya and C IV emission are C IV absorption line doublets. 

For most purposes normalized (continuum fitted) spectra from each order were used. 
The continuum fitting is made on a ±2000 km s" 1 region of interest by iterating a second- 
order polynomial fit to a 10-point smoothed spectrum and rejecting significantly deviant 
regions (|A(spectrum)/fit| > 0.025). The fit to the region around a C IV complex at z — 
3.26101 in Q1159+123 is shown in Figure Multiple Voigt profiles were then fitted to the 
absorption line profiles using the IDL version of Bevington's CURVEFIT algorithm. 
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3. Line Ratio Evolution with Redshift 
3.1. Introduction 

Each quasar spectrum was first searched for all complexes of C IV doublets redward of 
the La forest and outside of the quasar's proximity zone, defined to stretch to —4000 km s _1 
blueward of the quasar's redshift. The sample was then restricted to those cases for which 
SI IV also lies redward of the quasar's Lya emission. The redshift range used for each 
quasar is summarised in Table |l| Voigt profiles were next fitted to each such complex. The 
choice of the number of Voigt components is somewhat arbitrary because of the problems 
of blending but the complexes are quite weak in general and in the absence of saturation 
the integrated column densities of complexes should be well determined. In the same way, 
integrated column densities or upper limits for Si IV, C II and Si II were determined for 
lines lying longward of the forest. The usable redshift range for C II is also shown in 
Table |l|. The profile fits to the C IV and Si IV doublets in all the cloud complexes with 
10 13 cm~ 2 < iV(C IV) < 10 14 cm -2 are shown in Figure 5. The redshifts, number of Voigt 
components (a minimum of two was used), C IV column density, Si IV column density, and 
C II and Si II column densities (where these lines are redward of the forest) are summarised 
for all complexes with N(C IV) > 5 x 10 12 cm -2 in Table §. 



3.2. Evolution of total column density ratios 

Figure |6| shows the evolution of the integrated value of Si IV/C IV in the 76 complexes 
with 5 x 10 12 cm -2 < N(C IV) < 10 14 cm -2 as well as in the the higher redshift systems 
satisfying this condition in Savaglio et al.'s (1996) observations of Q0000— 263. As SC and 
Savaglio et al. (1996) previously noted, Si IV/C IV increases considerably between z = 2 
and z — 4, with a strong suggestion of a break in the distribution of values around z — 3. 
In contrast, the C II/C IV ratio shown in Figure [7] does not exhibit this behaviour: in 
most of the complexes where C II/C IV can be measured, it is weak (C II/C IV < 0.1) or 
undetected, regardless of redshift. Median values of Si IV/C IV are shown as a function of 
redshift in Figure |. The median ratio is 0.039±g;gg| for all the systems below z = 3 but 
0.13^o;o4 above z = 3, where the errors are 1 a computed using the median sign method. 
A key diagnostic in interpreting the data is whether the evolution is in the form of a step 
function or is a smoother change with redshift. A discontinuous change might arise from the 
effects of a change in the metagalactic ionizing spectrum at the He + ionization edge, whereas 
a smoother evolution might follow from the more general evolution of the external ionizing 
field and the range of internal densities and internal ionization in the complexes. The data 
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of Figure || and Figure || are strongly suggestive of an abrupt transition between z = 2.9 and 
z = 3. This redshift is slightly higher than, but comparable to, the redshift (z = 2.8 to 2.9) 
below which He II would be inferred to be predominantly ionized to He III from the direct 
observations of the He II Lya forest (Reimers et al. 1997). 

To quantify this, distribution functions of Si IV/C IV have been formed for complexes 
with z = 2.75 — 2.95, just below the suspected break, and for z = 3.0 — 3.2, just above 
it. The histograms of Si IV/C IV values for the 16 systems in the low redshift interval and 
the 22 systems in the high redshift interval are shown in Figure |9|. The median value of 
0.038 

-omo a t 2.75 < z < 2.95 is consistent with that of all the z < 3 systems, whereas the 
3.0 < z < 3.2 median value of 0.141q'o75 is consistent with that of all the z > 3 systems. 
Moreover, a rank sum test shows that the z = 2.75 — 2.95 data are consistent with being 
drawn from the same distribution as the z < 2.75 systems, and the z = 3.0 — 3.2 data with 
being drawn from the same distribution as the z > 3.2 systems. However, the rank sum test 
shows that the z = 3.0 — 3.2 data has a probability of only 2.2 x 10 -4 of being drawn from 
the same distribution as the z = 2.75 — 2.95 data, and that at the 98% confidence level the 
data require a jump of at least a multiplicative factor of 1.8 from the lower to the higher 
redshift interval. The data are therefore consistent with invariant distributions at high and 
low redshift together with an abrupt jump between z = 2.9 and z = 3. The corresponding 
step function, with Si IV/C IV changing by a factor of 3.4 at z — 2.95, is shown as the 
dashed line in Figure |[ 

Assuming instead that the evolution of the median ionization ratio is a smooth power 
law in (1 + z) requires an evolution of the form (1 + z) 5A which is shown as the dotted line in 
Figure §. This however predicts a jump of only 1.57 between z = 2.85 and z — 3.1 which, as 
already noted, is not large enough to reconcile the relative distribution functions of Figure [|. 
It would therefore seem that simple smooth power law evolution is not consistent with the 
data. 



3.3. Internal structure 

Variations of the ionization ratios in different regions within a complex can result from 
structure in the internal density of material in the complex, or from the presence of internal 
ionization sources. It is therefore of considerable interest to examine the ionization ratios 
within the complexes to understand how these effects can contribute to the overall ionization 
balance and to determine whether or not this can change the conclusions of section 3.2. 

The simplest approach to this is to break down the comlexes into the components 
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determined by Voigt profile fitting (Boksenberg 1998). Figure 10 shows the redshift evolution 
of the Si IV/C IV components broken down in this way. All components with iV(C IV) > 
5 x 10 12 cm -2 have been included irrespective of the total column density of the complex, 
so that a small number of high density complexes provide a large number of components, 
and these also include a number of Lyman limit systems in which radiative transfer effects 
are significant. However, even with these effects and with the increased confusion owing to 
the large spread in components, the change between high and low redshifts is clearly seen. 
Figure |ll| shows histograms of the distributions above and below z = 3. If the components are 
considered to be statisticaly independent samples then the median ratio is 0.043io]ool below 
z = 3 and 0.15^o;q4 above z = 3. At the 98% confidence level, the minimum multiplicative 
increase to bring the distributions into consistency is 1.93. 

However, the methodology of fitting Voigt profiles is a very unsatisfactory approach to 
this problem because of the considerable arbitrariness in the choice of components, which can 
have several unfortunate effects, including forcing other ions preferentially into or out of a 
particular component, blending components, and introducing non-independent components. 
These problems can be avoided by using the more direct technique of computing the ratios 
of the optical depths of various ions throughout the complex. This procedure is illustrated 
in Figures ^ and [13|. For each identified C IV complex, all data points lying above some 
limiting optical depth (here r(C IV 1550) > 0.05) were identified, and optical depths in the 
corresponding regions of the other ions were then measured. Si IV and C IV optical depths 



for a complex at z = 3.33714 in 2000—330 are shown in Figure [12[ The ratio (tA/) relative 
to that of C IV then gives the ratio of iV(ion)/iV(C IV) throughout the velocity range of the 
complex, as is illustrated in Figure [13], where it is compared with the average ion ratios in 
the complex, shown as the solid line. This procedure can be slightly in error if there is a 
substantial amount of thermal broadening but should be adequate for these heavier elements 
where the broadening is primarily kinematic. Si IV/C IV values determined in this way for 
all the complexes of Table ^| is shown in Figure [14] where again a clear evolution of the 
distribution can be seen. 

In order to determine the noise levels associated with this procedure, a large number of 
blank field realizations were constructed, in which the redshift was slightly displaced from 
the true value. This Monte Carlo technique is extremely powerful, providing an unbiased 
measure of the noise, including the systematic effects of continuum uncertainty, varying S/N 
among the quasars and within a single quasar spectrum, and contamination by unidentified 
absorption lines. Indeed, this last-named property allows a robust measure (§4) of the 
average column densities of ions lying in the Lyman alpha forest. 

The measured distribution (solid line) of Si IV/C IV is compared with the average of 
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the random realizations (dashed line) in Figure 15 for all the systems, and also for those 
above and below z = 3. Si IV/C IV is clearly seen to be much stronger above z = 3. 
Figure [T^ shows median Si IV/C IV values as a function of redshift, with errors derived from 
the distribution of medians in the blank-field realizations. The data again require a rapid 
jump at around z = 3, with values of 0.035 ± 0.005 below z = 3 and 0.10 ± 0.01 above. 



4. Ionization Ratios 
4.1. Introduction 

The local shape of the ionizing flux and the value of the ionization parameter may be 
estimated from the relative distribution of the ionization levels of the various species. Most 
work to date has concentrated on using Si IV/C IV versus C II/C IV, which, as we have 
discussed in the introduction, is a powerful discriminant of the ionization parameter and 
of the strength of any break near the He + edge. This diagnostic also relies on absorption 
lines that can be measured outside the forest. The Si IV/C IV versus C II/C IV method is 
reconsidered (§ 4.2) in the light of the present larger data set. However, the optical depth 
ratio method developed in §3.3 provides a powerful way to measure the average abundances 
of a number of other species (Si III, C III, N V, and O VI) that lie primarily within the Lya 
forest. We use this technique in §4.3 to determine the ensemble properties of these ions as 
a function of redshift and to show that they can be fitted into a consistent interpretation of 
the evolution of the shape of the metagalactic flux. 



4.2. Si IV/C IV versus C II/C IV 

Si IV/C IV versus C II/ C IV is shown in Figure |P7| for the 29 complexes with 10 13 cm -2 < 



N(C IV) < 10 14 cm -2 and for which the C II line lies longward of the Lya forest. Open 
squares denote systems with z > 3 and solid squares those with z < 3, while downward 
and/or leftward pointing arrows correspond to systems in which Si IV or C II is not detected, 
with the point positioned at the 1 a level. The curves show a range of models computed with 
the CLOUDY code (Ferland 1993) for a power-law ionizing spectrum with a range of breaks 
at the He + edge at 4 Rydberg, and with log 10 (Si/C) = —0.66, typical of metal-poor systems. 
As was found previously in SC, while many of the low redshift systems are consistent with 
photoionization by a pure power law, most of the high redshift systems have Si IV/C IV 
values which are much too high relative to C II/C IV for this to be the case. They are 
consistent with photoionization by heavily broken power laws with breaks across the He + 
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edge in excess of 20 and softness ratios, S = f v {& Ryd)// I/ (1 Ryd) in excess of 300. 

The models of SC were computed with simple broken power laws and a more extensive 
discussion of more general spectra may be found in Giroux & Shull (1997). However, Fig- 
ure [18| illustrates that the general behaviour is fairly model invariant; it compares the data 



with a wide range of power laws (—1.5 to —2.5) again computed with log 10 (Si/C) = —0.66. 
Unless Si/C is made unreasonably high at z > 3, not even the softest power law spectrum 
(the top curve) can approach the high Si IV/C IV values. The introduction of the break 
at 4 Ryd moves the curves to a new asymptote, as is illustrated by the dashed and dot- 
ted curves. Both are computed with an ionization parameter of —1.6 and a —2 power law 
spectrum below 4 Ryd, but with a range of break strengths. The dotted curve corresponds 
to the same models as in Figure [T7| in which the spectrum above the break is extrapolated 
with the same spectral index as it has below the break but at the lower flux level, as might 
be the case if there was a substantial galaxy contribution to the ionizing spectrum. The 
dashed line is computed instead with a model in which the break is a He + opacity whose 
cross section reduces as E~ 2 ' 7 at higher energies. This latter model asymptotes to the same 
Si IV/C IV values but requires a larger softness ratio of greater than about 2000 to achieve 
this. It is this effect that gave rise to Rauch et al.'s (1997) incorrect conclusion that He + 
breaks could not account for the high values of Si IV/C IV in high ionization systems since 
they used a recovering spectrum and an intermediate softness parameter (a flux reduction of 
100 at 4 Ryd). However, as the break becomes large, both models reach the same limiting 
values which are consistent with the observed high values of Si IV/C IV and low values of 
C II/C IV. 

In this asymptotic limit of a large break at the He + edge, the Si IV/C IV to C II/C IV 
ratio approaches a single track characterized essentially uniquely by the ionization parameter, 
which determines, modulo a very weak dependence on spectral index, the C II/C IV value. 
Therefore most of both the high and low redshift systems are consistent with a relatively small 
range of ionization parameter, log 10 T = — 1 to —2, but the higher redshift systems mostly 
require the high energy spectrum to be blocked above 54 eV whereas at the low redshifts 
this is not generally the case. It is perhaps important to emphasise that this effect is not 
observationally subtle and in Figure 19 we show C IV, Si IV and C II profiles for a sample of 
the high Si IV/C IV systems over the redshift range. For the high Si IV/C IV systems at low 
redshift, C II is generally clearly present but many of the high-z high Si IV/C IV systems 
are simply missing C II absorption. 
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4.3. Other ion ratios 

The optical depth technique discussed in section 3 can also be used to determine the 
average properties of other important ions such as C III, Si III and O VI whose transitions 
lie below the wavelength of Lya and for which individual line strengths cannot be reliably 
determined because of forest contamination. 

As with Si IV, the value of t(X)\xF x /t(C IV) Ac rv/c iv was determined for the ion, X, 
at all positions at which significant C IV absorption was present (again, r(C IV 1550) > 0.05) 
and a distribution of this quantity was constructed. A large number of random realizations 
were then made in which the C IV redshift was slightly displaced, to determine the back- 
grounds and errors. Because of forest contamination, which acts as an effective noise source, 
the errors are larger than for species that lie outside the forest, such as Si IV. 

Figure 20 shows the distributions of C III/C IV, N V/C IV and O VI/C IV for all 
the available systems, and separately for z > 3 and z < 3 systems where the division is 
possible. The median and mean values and their 1 a errors for a range of ions are given 
in Tables |3| and |] respectively. As is expected, O VI is strongly detected, with a median 
iV(0 VI)/N(C IV) = 1.97 ± 0.28. C III, Si III and N V are more weakly detected, with a 
median iV(C III)/JV(C IV) = 0.75 ± 0.25. Neither Si II nor C II is significantly detected. 

The ionization patterns above and below z = 3 are summarised in Figure 21 which 
uses the more robust median values. The points show the quantities A c N(X)/A X N(C IV), 
where (A x /Ac) is the assumed abundance of the element X with respect to C, which is 
taken to be [O/C] = 0.5, [Si/C] = 0.4 and [N/C] = —0.7, characteristic of low metallicity 
systems. This quantity is plotted against the ionization potential of the ion. 

The primary difference between z < 3 and z > 3 remains the much higher abundance 
of Si IV at the higher redshift. The fractions of the high ions such as N V and O VI are 
marginally larger at z < 3 than at z > 3, but the effect is not large. The z < 3 ionization 
distribution can be reproduced quite well by a very simple model in which the ionizing 



spectrum is a —1.8 power law and T = —1.6; this is shown as the solid bars in Figure |21a . 
This highly oversimplified model is of course only a rough representation of the ensemble 
average of conditions in the IGM. However, for z > 3 we need a mix of conditions. As is 
discussed above, most of the systems must lie in the regions exposed to highly broken spectra 
in order to reproduce the lower ion ratios, particularly Si IV/C IV, but some fraction of the 
volume must still be ionized by spectra containing higher energy photons to reproduce the 
observed O VI, if this is assumed to be produced by photoionization. A model in which 70% 
of the systems are irradiated by a fully broken spectrum with a = —1.8 below the break and 
T = —1.5, and 30% of the volume is occupied by a pure a = —1.8 spectrum with T = —1.0 



provides a reasonable representation, as is show by the bars in Figure 



5. Conclusions 

I summarise the results of the paper by noting that, irrespective of the analysis method- 
ology adopted, there is a significant change in the ionization balance of forest metal lines 
which occurs just below a redshift of 3. At lower redshifts, the ionization balance in the 
forest lines is fully consistent with a pure power law ionization spectrum with an index of 
— 1.8 but at higher redshifts the high values of Si IV/C IV seen in most of the forest clouds 
despite generally low C II/ C IV values implies that the ionizing flux must be very soft, with 
a large break at the He + edge. The change occurs quite rapidly between z = 2.9 and z — 3, 
just above the redshift at which highly patchy He II Lya absorption is seen in the quasar 
HE 2347—4342 (Reimers et al. 1997). The simplest explanation seems to be that we are see- 
ing the redshift at which He II ionizes completely to He III as the He III Stromgren spheres 
overlap. 

I am grateful to the many people at the Keck telescopes who assisted with these obser- 
vations and most particularly to Steve Vogt whose HIRES spectrograph made them possible. 
I would also like to thank Len Cowie and Esther Hu for obtaining some of the observations 
on which this work is based, Dieter Reimers for computing positions and magnitudes for un- 
published bright quasars, and Sandra Savaglio for providing updated information on higher 
redshift Si IV/C IV values. The research was supported by the National Science Foundation 
under grant AST 96-17216. 
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Table 1. Observations 



Quasar 


Mag. 


Expo, 
(hrs) 


S/N 


Complete 


Zlow 

(Si IV) 


Zlow 

(C II) 


Zhigh 


0014+813 


16.5 


7.65 


120 


Yes 


2.81 


3.00 


3.31 


0256-000 


18.2 


5.33 


60 


Yes 


2.78 


2.97 


3.28 


0302-003 


17.8 


10.67 


90 


Yes 


2.72 


2.91 


3.21 


0636+680 


16.5 


6.08 


130 


Yes 


2.62 


2.80 


3.10 


0741+4741 


16.1 


4.0 


145 


No 


2.66 


2.84 


3.15 


0956+122 


17.5 


6.67 


90 


Yes 


2.74 


2.92 


3.23 


1159+123 


17.5 


8.0 


85 


Yes 


2.90 


3.09 


3.42 


1422+231 


16.5 


9.3 


330 


Yes 


3.01 


3.21 


3.54 


1623+269 


16.0 


3.33 


60 


Yes 


2.12 


2.27 


2.48 


1700+6416 


16.1 


8.33 


260 


Yes 


2.27 


2.43 


2.68 


2000-330 


19.0 


4.67 


60 


Near 


3.14 


3.36 


3.69 


2126-158 


17.3 


4.45 


130 


No 


2.72 


2.90 


3.21 


2347-4342 


16.3 


6.67 


130 


Yes 


2.37 


2.54 


2.82 
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Table 2. C IV, Si IV, C II and Si II Column Densities 



Quasar z # cpts. JV(C IV) JV(Si IV) 7V(C II) JV(Si II) 



1 623+269 


9 1 61 41 

Zj . lul^tl 


•3 

O 


6 1 

W . X 


y 
a 


in 13 

X w 


1 4 

X . t: 


v- 

A 


in 12 

X w 






— 1 9 

X . u 


y 

A 


in 11 

X w 


1 623+269 


9 24450 

Z .Z 1 1(JU 


2 

Zj 


3 1 

W . X 


x 


10 13 

X w 


—3 2 

O . Zj 


x 


10 10 

X w 






—2 4 

Zj .T: 


y 

A 


10 11 

X W 


1 7(10+641 6 


9 98Q59 

Zj . ZjO<_/ WZj 


9 

Zj 


6 8 
w .o 


v 
a 


10 12 

X w 


6 9 

U. Zj 


v 

A 


10 10 

X w 






—9 6 

Zj .U 


y 

A 


10 11 

X w 


1700+6416 


2 30826 


W 


3 4 

W . T 


x 


10 13 

X w 


1 7 

X . 1 


x 


10 12 

X w 






—9 5 

U • 'J 


x 


10 11 

X w 


X i WW \^\Jrt X W 


2 31 536 

Zj . O X WOW 


1 3 


1 9 

X . U 


V 
A 


10 15 

x w 


6 1 

W. X 


A 


10 13 

X w 






9 5 

U . O 


y 

A 


10 12 

X w 


x i wwh^wtixw 


2 37984 

Zj . W 1 l/Ut: 


2 

Zj 


1 3 

X . O 


y 
a 


10 13 

X w 


7 3 

I . W 


y 

A 


10 11 

X w 






4 6 

T . W 


y 

A 


10 11 

X w 


1 623+269 


2 40085 

Zj .tiWWwW 


5 

W 


1 2 

X . Zj 


v- 

A 


10 14 

X w 


— 1 8 

X . O 


y 

A 


10 11 

X w 


9 7 x 


10 10 

X w 


4 Fi 

t: . W 


y 

A 


10 11 

X w 


1 7nf)-i-fi41 6 


9 43304 


7 


4 3 

T . - J 


A 


10 13 

X w 


Q (1 

C/. W 


y 

A 


10 12 

X w 






— 1 9 

X . Zj 


y 

A 


10 11 

X w 


9347-4349 


9 43814 

Zj .4:00 It: 


9 

Zj 


7 5 

1 . W 


v 

A 


10 12 

X w 


3 (1 

O. W 


y 

A 


10 11 

X w 






—3 4 

O . t 


y 

A 


10 10 

X w 


1 7(1(1-1-641 6 

x i wwh^wtixw 


2 43863 


fi 

O 


5 2 

W . Zj 


x 


10 13 

X w 


2 

Zj . W 


x 


10 12 

X w 


2 6 x 

Zj . W A 


10 11 

X w 


— 1 3 

X . O 


x 


10 12 

X w 


1 623+269 

IUZjO |^ Zj W t7 


2 44353 


W 


1 7 

X . 1 


x 


10 13 

X w 


6 8 

W . O 


x 


10 10 

X w 


— 3 1 x 

W . X A 


10 12 

X w 


—1 2 

X . Zj 


x 


10 12 

X w 


1 693+96Q 


9 44534 

Z, . tit: W -Jt: 


9 

Zj 


6 3 

W . W 


v 
a 


10 12 

X w 


1 (1 

X . w 


y 

A 


10 12 

X w 


9 7 x 

Zj . 1 A 


10 12 

X w 


1 4 

X . t: 


y 

A 


10 10 

X w 


1 7(1(1-1-641 6 


9 5681 7 

Zv . WWO X l 


9 

Zj 


6 6 

W . W 


v 

A 


10 12 

X w 


9 6 

Zj. U 


y 

A 


10 11 

X w 


— 5 5 x 

O.O A 


10 10 

X w 


9 (1 

Zj. W 


y 

A 


10 10 

X w 


1 7(1(1-1-641 6 

X 1 WWh^UtiXU 


9 57857 

Z.tJ 1 OW 1 


4 

t: 


3 6 

W . W 


v 

A 


10 13 

X w 


1 3 

X . O 


y 

A 


10 12 

X w 


1 Q x 

X . <_/ A 


10 12 

X w 


—3 6 

O . W 


y 

A 


10 11 

X w 


2347-4342 

Z. W t: ( t: -.J t: Z. 


2 63449 

Zj ,UU 1 1J 


9 

Z 


9 9 


A 


10 12 

X w 


1 8 

X . O 


y 

A 


10 10 

X w 


Q Q x 

U . U A 


10 12 

X w 


—3 

o . w 


y 

A 


10 11 

X w 


(1741 -1-4741 


2 67283 

Zj . W 1 zoo 


2 

Zj 


1 

X . w 


v- 

A 


10 13 

X w 


8 7 


y 

A 


10 10 

X w 






1 3 

X . o 


y 

A 


10 11 

X w 


0636+68(1 


9 68998 

Zj .uozzo 


9 

Zj 


8 4 




10 12 

X w 


7 8 


y 

A 


10 11 

X w 






1 (1 

X . w 


y 

A 


10 12 

X w 


(1741 -1-4741 

W 1 t: x n^^t ' x 


9 6Q935 

Zj .UC/ZjWW 


9 

Zj 


6 5 

U . O 


A 


10 13 

X w 


4 (1 

L ±. W 


y 

A 


10 12 

X w 






—9 7 

Zj . 1 


y 

A 


10 11 

X w 


(1741 -1-4741 


2 69478 

Zj . W Jt: 1 (J 


9 

Zj 


2 7 

Zj . 1 


x 


10 13 

X w 


6 8 

W . O 


x 


10 11 

X w 






— F5 Fi 

O . O 


x 


10 11 

X w 


(1741 -1-4741 


2 71446 

Zj . I X J. J.W 


9 

Z 


6 3 

W . O 


x 


10 12 

X w 


2 8 

Zj . (J 


x 


10 11 

X w 






1 9 

X . U 


y 

A 


10 11 

X w 


91 96-1 58 


2 72796 




2 8 


x 


10 13 

X w 


2 3 


x 


10 12 

X w 






—6 1 


y 

A 


10 11 

X w 


U ( 41+4 ( 41 


O 7QQOQ 


Z 


0.1 


X 


1 n!3 

1U 


O.O 


X 


1 n 12 

1U 






1 1 
— 1.1 


X 


m 12 

1U 


2347-4342 


2.73564 


5 


1.1 


X 


10 14 


1.9 


X 


10 13 


6.8 x 


10 12 


4.8 


X 


10 10 


2126-158 


2.76916 


12 


4.6 


X 


10 14 


6.4 


X 


10 13 






6.5 


X 


10 14 


0741+4741 


2.78488 


2 


1.0 


X 


10 13 


3.9 


X 


10 11 






-1.6 


X 


10 11 


0302-003 


2.78663 


9 


7.1 


X 


10 13 


3.1 


X 


10 12 






4.9 


X 


10 11 


2126-158 


2.81959 


5 


3.0 


X 


10 13 


2.0 


X 


10 12 






1.7 


X 


10 11 


0256-000 


2.82444 


6 


2.3 


X 


10 13 


1.2 


X 


10 12 






-7.1 


X 


10 11 


0302-003 


2.82593 


2 


6.8 


X 


10 12 


-2.4 


X 


10 11 






2.1 


X 


10 11 


0956+122 


2.83185 


2 


1.3 


X 


10 13 


8.3 


X 


10 10 






-6.9 


X 


10 11 


0956+122 


2.83433 


2 


8.9 


X 


10 12 


2.3 


X 


10 11 






-2.1 


X 


10 11 


0256-000 


2.83599 


2 


5.6 


X 


10 12 


4.2 


X 


10 11 






8.0 


X 


10 11 


0636+680 


2.86883 


2 


7.3 


X 


10 12 


5.0 


X 


10 11 


-3.9 x 


10 11 


4.7 


X 


10 11 


0956+122 


2.88806 


4 


1.3 


X 


10 13 


5.1 


X 


10 11 






1.1 


X 


10 12 
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Table 2 — Continued 



Quasar z # cpts. N(C IV) JV(Si IV) N(C II) JV(Si II) 



UUUU | uuu 


2 89163 


5 


2.8 


x 


10 13 


5.6 


x 


10 10 


4.5 x 10 9 


-6.6 x 10 11 


0636+680 


2 90374 


9 


3.0 


x 


10 14 


1.4 


x 


10 14 


4.3 


x 


10 14 


5.5 x 10 13 


0741+4741 


2 90455 


6 


2.4 


x 


10 13 


7.8 


x 


10 11 


2.2 


x 


10 12 


—1.4 x 10 12 


2126-158 


2 90708 

• kJ \J 1 WW 


2 


3.6 


x 


10 13 


8.7 


x 


10 12 


5.5 


X 


10 12 


-3.1 x 10 11 


0014+813 

W W X j: i X W 


2 90823 


2 


6.8 


x 


10 12 


4.5 


x 


10 11 








6.5 x 10 11 


0956+122 

W W W \ X ^1 ^1 


2 91490 


2 


8.7 


x 


10 12 


2.1 


x 


10 11 








-3.4 x 10 s 


0302-003 


2.91749 


2 


1.8 


x 


10 13 


5.0 


x 


10 11 


1.3 


X 


10 12 


4.0 x 10 11 


0302-003 

VJUUZj WWW 


2 95882 


2 


7.2 


x 


10 12 


1.1 


x 


10 11 


-1.1 


x 


10 12 


—3.5 x 10 11 


2126-158 

i—t X i—t \J X W^J 


2 96346 

• • J V.' lVI 


3 


2.0 


x 


10 13 


5.5 


x 


10 12 


1.6 


x 


10 11 


9.8 x 10 11 


0741+4741 


2 96530 

• i_/ W W W W 


3 


1.0 


x 


10 13 


—1.3 


x 


10 11 


2.6 


x 


10 11 


6.7 x 10 11 


2126-158 


2.96747 


5 


2.5 


x 


10 13 


3.8 


x 


10 12 


9.1 


X 


10 12 


7.0 x 10 11 


0302-003 

WWW^J WWW 


2 99540 


2 


6.6 


x 


10 12 


3.2 


x 


10 11 


-1.9 


X 


10 12 


-2.4 x 10 11 


0302-003 

\JKJ\J i-t WWW 


3 00295 


2 


3.3 


x 


10 13 


1.7 


x 


10 12 


7.4 


x 


10 11 


1.2 x 10 12 


0Q56+122 


3 01043 


2 


9.7 


x 


10 12 


1.5 


x 


10 12 


3.8 


X 


10 12 


—1.0 x 10 12 


0636+680 

X/ W W W | VJU V 


3 01303 


2 


1.3 


x 


10 13 


5.8 


x 


10 11 


-1.4 


X 


10 12 


-7.8 x 10 11 


0636+680 

WWWW \ VJU V 


3.01749 


2 


4.1 


x 


10 13 


1.2 


x 


10 12 


7.0 


X 


10 11 


-4.1 x 10 9 


0741+4741 


3 01763 

W • W X 1 V 7 7 


6 


6.6 


x 


10 13 


2.5 


x 


10 13 


1.1 


X 


10 15 


1.5 x 10 14 


0256-000 

WZvWW WWW 


3 01799 


2 


1.5 


x 


10 13 


1.7 


x 


10 12 


4.3 


x 


10 10 


—1.2 x 10 12 


0741+4741 

W I j: X | i 1 l. X 


3 03469 


3 


1.0 


x 


10 13 


9.7 


x 


10 11 


5.9 


x 


10 10 


—2.3 x 10 11 


0302-003 

uUUi/ WWW 


3 04705 

"7 . V 7 l 1 WW 


2 


5.7 


x 


10 12 


-8.7 


x 


10 10 


3.2 


X 


10 11 


—5.7 x 10 10 


0956+122 

\j <^> w w \ x — - — - 


3 05287 


2 


1.0 


x 


10 13 


1.8 


x 


10 12 


-3.3 


X 


10 11 


6.1 x 10 10 


0741+4741 


3 05366 

W • W l/U w w 


4 


1.1 


x 


10 13 


6.3 


x 


10 12 


1 7 


A 


1U 


1.7 x 10 11 


UZOO-UUU 


q nfi/i Q/i 


o 

z 


4.0 


X 


1 n!3 


1 Q 
l.O 


X 


1 n!3 

1U 


6.6 


X 


10 12 


i i \s i nil 

L.I X 1U 


1422+231 


3.09030 


4 


4.2 


X 


10 13 


3.1 


X 


10 12 








-1.3 x 10 11 


0956+122 


3.09691 


2 


7.1 


X 


10 12 


1.0 


X 


10 12 


1.5 


X 


10 12 


8.6 x 10 11 


0956+122 


3.11418 


2 


6.9 


X 


10 13 


3.9 


X 


10 12 


1.2 


X 


10 12 


-5.2 x 10 11 


1422+231 


3.13413 


3 


1.6 


X 


10 13 


9.7 


X 


10 11 








-3.3 x 10 10 


1422+231 


3.13702 


2 


6.2 


X 


10 12 


1.1 


X 


10 12 








5.3 x 10 11 


0956+122 


3.15294 


5 


3.9 


X 


10 13 


1.2 


X 


10 13 


4.2 


X 


10 12 


3.3 x 10 11 


1159+123 


3.16701 


2 


8.6 


X 


10 12 


5.5 


X 


10 11 


-3.6 


X 


10 11 


1.5 x 10 11 


2000-330 


3.17249 


2 


1.3 


X 


10 13 


6.7 


X 


10 12 








2.1 x 10 12 


0956+122 


3.17860 


3 


1.6 


X 


10 13 


2.3 


X 


10 12 


8.8 


X 


10 11 


-3.0 x 10 11 


2000-330 


3.19155 


5 


7.6 


X 


10 13 


1.9 


X 


10 13 








6.2 x 10 13 


0256-000 


3.19849 


2 


1.6 


X 


10 13 


2.4 


X 


10 11 


6.3 


X 


10 11 


2.6 x 10 11 



Table 2 — Continued 



Quasar z # cpts. JV(C IV) JV(Si IV) JV(C II) JV(Si II) 



0956+122 


3.22310 


5 


7.5 


X 


10 13 


1.5 


X 


10 13 


2.0 


X 


10 12 


1.9 x 10 11 


1159+123 


3.22554 


2 


2.6 


X 


10 13 


2.2 


X 


10 12 


-1.3 


X 


10 12 


-2.7 x 10 11 


0014+813 


3.22623 


5 


5.8 


X 


10 13 


3.0 


X 


10 12 


-1.6 


X 


10 12 


-1.3 x 10 12 


0014+813 


3.23339 


2 


9.6 


X 


10 12 


1.0 


X 


10 12 


2.0 


X 


10 11 


7.0 x 10 10 


1159+123 


3.26101 


2 


3.8 


X 


10 13 


2.6 


X 


10 12 


-4.0 


X 


10 11 


7.6 x 10 11 


2000-330 


3.33281 


4 


6.2 


X 


10 13 


2.9 


X 


10 13 








3.8 x 10 12 


2000-330 


3.33714 


2 


2.7 


X 


10 13 


3.7 


X 


10 12 








8.2 x 10 11 


1159+123 


3.37844 


2 


1.6 


X 


10 13 


2.2 


X 


10 12 


1.3 


X 


10 11 


-5.2 x 10 9 


1422+231 


3.38277 


5 


4.4 


X 


10 13 


2.8 


X 


10 12 


2.1 


X 


10 12 


3.7 x 10 11 


1159+123 


3.40448 


2 


1.9 


X 


10 13 


2.3 


X 


10 12 


-6.8 


X 


10 11 


3.5 x 10 11 


1422+231 


3.41121 


2 


8.7 


X 


10 12 


1.7 


X 


10 12 


4.7 


X 


10 10 


7.9 x 10 10 


1422+231 


3.44717 


2 


4.0 


X 


10 13 


7.6 


X 


10 12 


-4.0 


X 


10 10 


3.9 x 10 10 


2000-330 


3.50532 


2 


9.7 


X 


10 12 


4.1 


X 


10 12 


2.6 


X 


10 12 


2.3 x 10 10 


1422+231 


3.53863 


11 


1.8 


X 


10 14 


3.8 


X 


10 13 


4.4 


X 


10 13 


2.4 x 10 12 


2000-330 


3.54937 


11 


2.3 


X 


10 14 


6.3 


X 


10 13 


1.7 


X 


10 14 


9.4 x 10 12 
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Table 3. Medians of Column Density Ratios With Respect to C IV 



Ion 


All 


z > 3 


z < 3 


C II 


0.02 ±0.03 


0.01 ±0.04 


0.03 ±0.09 


C III 




0.75 ±0.25 




C IV 


1 


1 


1 


Si II 


-0.001 ±0.004 


0.004 ±0.010 


-0.004 ±0.006 


Si III 


0.029 ±0.010 


0.08 ±0.037 


0.013 ±0.019 


Si IV 


0.056 ±0.003 


0.10 ±0.005 


0.035 ±0.004 


N V 


0.019 ±0.008 


0.018 ±0.024 


0.019 ±0.006 


VI 


1.97 ±0.28 


1.59 ±0.21 


2.75 ± 0.78 
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Table 4. Means of Column Density Ratios With Respect to C IV 



Ion 


All 


z > 3 


z < 3 


C II 


-0.28 


-0.03 


-0.47 


C III 




1.97 




C IV 


1 


1 


1 


Si II 


0.08 


0.27 


-0.07 


Si III 


0.0 


0.07 


-0.04 


Si IV 


0.12 


0.17 


0.07 


N V 


0.05 


0.09 


0.02 


VI 


1.8 


1.3 


2.2 




5620 5640 5660 5680 5700 5720 

VACUUM HELIOCENTRIC WAVELENGTH (A) 



Fig. la. — A sample extraction of one of the orders longward of the Lya forest in the quasar Q1159+123. The total 
exposure is 2 hr, consisting of three 40 m integrations. The spectrum is unsmoothed and oversamples the resolution 
element. 
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Fig. lb. — An order in the same spectrum as (a) but lying in the Lja forest. The precision of the sky subtraction may 
be assessed from the degree to which the bases of the saturated Lya lines approximate to zero. 
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Fig. 2. — Final processed and extracted order in Q1159+123 for comparison with Figure 1(a). The total exposure of 
8 hr is combined from a number of sub-exposures with spectra covering each order added in a S/N optimized fashion. 
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Fig. 3. — Complete spectrum of Q1159+123 formed by interpolating all the extracted orders to a uniform 0.05 A 
wavelength scale and then splicing adjacent orders. The principal feature is the quasar's Lya emission near 5400 A and 
the spectrum extends just past C IV emission near 7000 A. There is a Lyman limit system near 4150 A. Most of the 
features longward of the Lya emission are C IV doublets, except for the atmospheric A band. 
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Fig. 4. — A typical continuum-fitted region of spectrum, around the z = 3.26101 C IV absorption line complex 
in Q1159+123. The continuum was fitted to a ±2000 km s" 1 region around the C IV A1548 A line by iterating a 
2nd-order polynomial fit to a 10-pt smoothed spectrum and rejecting significantly deviant regions. 
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Fig. 5a-i. — Voigt profile fits to the C IV and Si IV doublets in all cloud complexes in the sample with 10 13 cm~ 2 
N(C IV) < 10 14 cm~ 2 . All doublets lie redward of the quasar's Lya emission. Table 2 summarises the redshifts, numb 
of Voigt components, total C IV and Si IV column densities of each complex, and C II and Si II column densities. 
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Fig. 6. — The ratio of total Si IV column density to total C IV column density as a function of redshift for the 76 
complexes with 5 x 10 12 cm" 2 < N(C IV) < 10 14 cm" 2 (filled squares). Open squares denote values from Savaglio et 
al.'s (1997) observations of Q0000— 263. Downward pointing arrows show systems in which Si IV is not detected and 
where the points are positioned at the 1 a upper limit. 
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Fig. 7. — The ratio of total C II column density to total C IV column density as a function of redshift for the complexes 
of Figure 6 in which the C II system lies longer in wavelength than the forest. Downward pointing arrows show systems 
in which C II is not detected and where the points are positioned at the 1 a upper limit. 
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Fig. 8. — Median values of Si IV/C IV as a function of redshift for the sample of Figure 6. Errors are ±1 o computed 
using the median sign method. The median value of Si IV/C IV is 0.039±g;gg§ for z < 3 and 0.13lo;o4 f° r z > 3. The 
dashed line is a step function fit to the data, with a jump of a factor of 3.4 at z — 2.95. The dotted line has the form 
(1 + zf A . It does not predict a sufficient jump between z = 2.85 and z — 3.1 to reconcile the low and high redshift 
distributions (Figure 9). 
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Fig. 9. — Histograms of the distributions of log 10 (Si IV/C IV) for 22 systems with z = 3.0 — 3.2 (top) with a median 
of 0.14! t Ig;o7 5 and 16 systems with z = 2.75 — 2.95 (bottom) with a median of 0.038!t(xoio. A multiplicative jump of a 
least a factor of 1.8 is required at the 98% confidence level to reconcile the low and high redshift distributions. 
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Fig. 10. — Si IV/C IV as a function of redshift for all Voigt profile components with iV(C IV) > 5 x 10 12 cm 2 in the 
sample, irrespective of the total column density of the complex, and including a number of Lyman limit systems. 
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Fig. 11. — Histograms of the distributions of log 10 (Si IV/C IV) for all Voigt profile components in the sample of 
Figure 10 with z > 3 (top) with a median of 0.15lo!o4 an d for z < 3 (bottom) with a median of 0.043lo;oos- A 
multiplicative jump of a least a factor of 1.93 is required at the 98% confidence level to reconcile the low and high 
redshift distributions. 




Fig. 12.— Sample optical depths of Si IV (top) and C IV (bottom) for the z = 3.33714 complex in Q2000-330. The 
Si IV optical depth has been offset by 1.5. 
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Fig. 13. — Si IV/C IV values across the z = 3.33714 complex in Q2000— 330, computed from the Si IV and C IV optical 
depths of Figure 12 with r(C IV) > 0.05, as iV(ion) = r\f, where A is the wavelength and / the oscillator strength 
of the transition. The solid line is the average value of Si IV/C IV in the complex, computed from the total column 
densities of Si IV and C IV measured from Voigt profile fitting. 
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Fig. 14. — The Si IV/C IV density ratios computed from optical depths as in Figure 13 for all the absorption line 
complexes of Table 2, as a function of the redshift of the complex. 
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Fig. 15a-c. — The distribution of Si IV/C IV density ratios for all the data of Figure 14 (solid line) compared to the 
average of random realizations (dashed line). Median and mean values and 1 a errors are given in Tables 3 and 4, 
respectively. 
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Fig. 15b. — As in (a) for z > 3. 
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Fig. 15c. — As in (a) for z < 3. 
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Fig. 16. — Median Si IV/C IV as a function of redshift for the data of Figure 14, with ±1 a errors. The median is 
0.035 ± 0.005 for z < 3 and 0.10 ± 0.01 for z > 3. 
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Fig. 17.— Si IV/C IV versus C II/C IV for complexes with 10 13 cm- 2 < N(C IV) < 10 14 cirr 2 and with C II longward 
of the forest for systems with z < 3 (filled squares) and z > 3 (open squares). Arrows mark systems in which Si IV or 
C II is not detected, with the square positioned at the 1 a level. The curves are models computed with the CLOUDY 
code (Ferland 1993) with log 10 (Si/C) = —0.66 for a —1.5 power law ionizing spectrum (dashed line) and with the 
addition of a break at the He + edge at 4 Ryd of a factor of 2 (dotted line), 10 (dash-dot line), 100 (long dash line) and 
1000 (solid line). 
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Fig. 18. — Si IV/C IV versus C II/C IV for the data of Figure 17 (filled and open squares). The solid lines are 
CLOUDY models with power law spectra with index —1.5 (bottom), —2.0 (middle), and —2.5 (top), computed with 
log 10 (Si/C) = —0.66. The dashed and dotted curves show the effect of the introduction of a break in the spectrum at 
the He + edge. See text (§4.2) for details. 
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Fig. 19a. — Sample C IV, Si IV and C II profiles for systems with a high value of Si IV/C IV at z < 3. 
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Fig. 19b. — As in (a) for systems with z > 3, illustrating the absence of C II in these high redshift systems. 
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Fig. 20a-g. — (Solid lines): distributions of C III/ C IV, N V/C IV, and O VI/C IV computed from the optical depth 
ratios for r(C IV) > 0.05, as in Figure 15, compared with average blank realizations (dashed lines). Where possible, 
distributions are also shown separately for z > 3 and z < 3. Median and mean values and 1 a errors are given in 
Tables 3 and 4, respectively. 
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Fig. 21a. — The quantity Ac N(X)/Ax N(C IV) for z < 3 for the ions, X, of Table 3, where (Ax/ Ac) is the abundance 
of the element X with respect to C, as a function of the ionization potential of the ion, in eV. Errors are ±1 a. The 
solid bars show a model in which the ionizing spectrum is a —1.8 power law with T = —1.6. 
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Fig. 21b. — As in (a) for z > 3. See text (§4.3) for a description of the comparison model (solid bars). 



